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I We briefly review the diagonalization of quantum Hamiltonians using the quasi- 

sparse eigenvector (QSE) method. We also introduce the technique of stochastic 
^ . error correction, which systematically removes the truncation error of the QSE result 

J^l^ I by stochastically sampling the contribution of the remaining basis states. 
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Q I 1 Quasi-sparse eigenvector method 

S I Quasi-sparse eigenvector (QSE) diagonalization is a new computational method 

I which finds approximate low-lying eigenvalues and eigenvectors for a general 

quantum field Hamiltonian H [1]. It handles the exponential increase in the 
CJ ■ dimension of Fock space by exploiting the sparsity of the Hamiltonian. The 

method is most effective when the splitting between low-lying eigenvalues is 
^ I not too small compared to the size of the Hamiltonian off-diagonal entries, 

i-j I In such cases the low-lying eigenvectors are quasi-sparse, meaning that the 

vector is dominated by a small fraction of its largest components. The QSE 
algorithm is applied using the following steps: 

(1) Select a subset of ortho normal basis vectors {cj^,--- ,ei^} and call the 
corresponding subspace S. 

(2) Diagonalize H restricted to S and find one eigenvector v. 

(3) Sort the basis components {ei.\v) according to their magnitude and re- 
move the least important basis vectors. 

(4) Replace the discarded basis vectors by new basis vectors. These are 
selected at random from a pool of candidate basis vectors which are 
connected to the old basis vectors through non- vanishing matrix elements 
of H. 

(5) Redefine S as the subspace spanned by the updated set of basis vectors 
and repeat steps 2 through 5. 
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If the subset of basis vectors is sufficiently large, the exact low energy eigen- 
vectors will be stable fixed points of the QSE update process. We can show 
this as follows. Let \i) be the eigenvectors of the submatrix of H restricted to 
the subspace S, where S is the span of the subset of basis vectors after step 3 
of the QSE algorithm. Let \Aj) be the remaining orthonormal basis vectors 
in the full space not contained in S. We can represent H as 
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{1\H\A^) {1\H\A2) 
{2\H\A^) {2\H\A2) 



{A,\H\1) {A,\H\2) ■•• E-Xa, {Ai\H\A2) 
(A2I H\l) {A2\H\2) ■■■ {A2\ H \A,) E ■ Xa, 



(1) 



We have written the diagonal terms for the basis vectors \Aj) with an explicit 
factor E for reasons to be explained shortly. We let |1) be the approximate 
eigenvector of interest and shift the diagonal entries so that Ai = 0. Our 

starting hypothesis is that |1) is close to some exact eigenvector of H which 
we denote as |lfuii). More precisely we assume that the components of |lfuii) 
outside S are small enough so that we can expand in inverse powers of the 
introduced parameter E. 



We now expand the eigenvector as 



LfuU/ 



^mE-' + 
c'a.E-' + 



(2) 



and the corresponding eigenvalue as 

Afuii = a;^;-^ + • • • . (3) 

In (2) we have constrained |lfuii) to have unit norm. Prom the eigenvalue 
equation H jlfuii) = Afuu |lfuii) we find at lowest order 
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We see that at lowest order the component of |lfuii) in the \Aj) direction is 
independent of the other vectors 1^^'). If |1) is sufficiently close to |lfuii) 
then the limitation that only a fixed number of new basis vectors is added 
in step 4 of the QSE algorithm is not important. At lowest order in 
the comparison of basis components in step 3 (in the next iteration) is the 
same as if we had included all remaining vectors \Aj) at once. Therefore at 
each update only the truly largest components are kept and the algorithm will 
converge to some optimal approximation of | Ifun) . This is consistent with the 
actual performance of the algorithm for a wide range of examples. 



2 Stochastic error correction 



The quasi-sparse eigenvector method approximates the low energy eigenvalues 
and eigenvectors of a Hamiltonian H by finding and diagonalizing optimized 
subspaces ^S". This procedure however has an intrinsic error due to the omission 
of basis states not in S. The goal of stochastic error correction is to sample the 
contribution of these remaining basis vectors in a systematic fashion. Due to 
space limitations we consider here only the series expansion method. Several 
other techniques are discussed in a forthcoming paper [2] . 

The series method for stochastic error correction is based on the expansion 
introduced in the previous section, 



|lfull) OC 



1 

c'^E-^ + 4E-^ + ■ 

c'a^E-' + c':,,E-' + 
c'.E-^ + c'LE-' + 



(5) 



Xi^n = KE'' + KE~''--- . (6) 

This time it is more convenient to choose the normalization of the eigenvector 

such that the |1) component remains 1. The convergence of the expansion 
is controlled by the proximity of |1) to |lfuii). If |1) is not at all close |lfuii) 
then it will be necessary to use a more robust method, and such methods are 
discussed in [2]. 
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At first order in £" we find 



Aa, 



(7) 



{l\H\A_i}{A_i\H\l) 



(8) 
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{j\H\Ak){Ak\H\l} 



(9) 



The second order contributions are 



-Z^ Aa,. Aaj. •i--' .^-^ 



{A,\H\l){l\H\A^,){A^:\H\l) 



(10) 
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{\\H\Ai){Ai\H\A^)^A^\H\\) 
Aa^Aaj^ 



EEE 



m\Ai){A^\H\l){l\H\A^){A^\H\\) 
AAjAiAAj, 



+ EEE 



(m|if|A,)(A,|if|Afc)(Afc|g|l) 
Aaj- Aaj^, 



{m\H\A^){A,\H\l)m\A^){Ak\H\\) 

AmAAj '^''^'4fe 



E 



(l|f/|Aj)(Aj|J/|l) 
Aa, 



E 

L fe 



(m|f/|Afc)(Afc|g|l) 

A At, 
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These contributions are calculated by a Monte Carlo samphng process. Let 
-P(^triai(j)) denote the probability of selecting Atriai(i) on the i**^ trial. For 
example the first order correction to the eigenvalue can be calculated using 



a; = -e 



(l|g|A,){A,|g|l) 
Aa, 



lim -i: 

1=1,-, N 



(ll-^^l^tria^i)) (Arial(i) \H\i) 
•^^trial(i)^(^fial(i)) 



(13) 
(14) 



Application of the series method for stochastic error correction to the Hubbard 
model was discussed by Nathan Salwen at the Heidelberg meeting and is 
presented in another article in these proceedings. 
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